\(\int \frac {x^2 \sin (c+d x)}{(a+b x)^3} \, dx\) [34]

Optimal result

Integrand size = 17, antiderivative size = 241 \[ \int \frac {x^2 \sin (c+d x)}{(a+b x)^3} \, dx=-\frac {a^2 d \cos (c+d x)}{2 b^4 (a+b x)}-\frac {2 a d \cos \left (c-\frac {a d}{b}\right ) \operatorname {CosIntegral}\left (\frac {a d}{b}+d x\right )}{b^4}+\frac {\operatorname {CosIntegral}\left (\frac {a d}{b}+d x\right ) \sin \left (c-\frac {a d}{b}\right )}{b^3}-\frac {a^2 d^2 \operatorname {CosIntegral}\left (\frac {a d}{b}+d x\right ) \sin \left (c-\frac {a d}{b}\right )}{2 b^5}-\frac {a^2 \sin (c+d x)}{2 b^3 (a+b x)^2}+\frac {2 a \sin (c+d x)}{b^3 (a+b x)}+\frac {\cos \left (c-\frac {a d}{b}\right ) \text {Si}\left (\frac {a d}{b}+d x\right )}{b^3}-\frac {a^2 d^2 \cos \left (c-\frac {a d}{b}\right ) \text {Si}\left (\frac {a d}{b}+d x\right )}{2 b^5}+\frac {2 a d \sin \left (c-\frac {a d}{b}\right ) \text {Si}\left (\frac {a d}{b}+d x\right )}{b^4} \]

[Out]

-2*a*d*Ci(a*d/b+d*x)*cos(-c+a*d/b)/b^4-1/2*a^2*d*cos(d*x+c)/b^4/(b*x+a)+cos(-c+a*d/b)*Si(a*d/b+d*x)/b^3-1/2*a^
2*d^2*cos(-c+a*d/b)*Si(a*d/b+d*x)/b^5-Ci(a*d/b+d*x)*sin(-c+a*d/b)/b^3+1/2*a^2*d^2*Ci(a*d/b+d*x)*sin(-c+a*d/b)/
b^5-2*a*d*Si(a*d/b+d*x)*sin(-c+a*d/b)/b^4-1/2*a^2*sin(d*x+c)/b^3/(b*x+a)^2+2*a*sin(d*x+c)/b^3/(b*x+a)

Rubi [A] (verified)

Time = 0.38 (sec) , antiderivative size = 241, normalized size of antiderivative = 1.00, number of steps used = 14, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.294, Rules used = {6874, 3378, 3384, 3380, 3383} \[ \int \frac {x^2 \sin (c+d x)}{(a+b x)^3} \, dx=-\frac {a^2 d^2 \sin \left (c-\frac {a d}{b}\right ) \operatorname {CosIntegral}\left (x d+\frac {a d}{b}\right )}{2 b^5}-\frac {a^2 d^2 \cos \left (c-\frac {a d}{b}\right ) \text {Si}\left (x d+\frac {a d}{b}\right )}{2 b^5}-\frac {a^2 d \cos (c+d x)}{2 b^4 (a+b x)}-\frac {a^2 \sin (c+d x)}{2 b^3 (a+b x)^2}-\frac {2 a d \cos \left (c-\frac {a d}{b}\right ) \operatorname {CosIntegral}\left (x d+\frac {a d}{b}\right )}{b^4}+\frac {2 a d \sin \left (c-\frac {a d}{b}\right ) \text {Si}\left (x d+\frac {a d}{b}\right )}{b^4}+\frac {\sin \left (c-\frac {a d}{b}\right ) \operatorname {CosIntegral}\left (x d+\frac {a d}{b}\right )}{b^3}+\frac {\cos \left (c-\frac {a d}{b}\right ) \text {Si}\left (x d+\frac {a d}{b}\right )}{b^3}+\frac {2 a \sin (c+d x)}{b^3 (a+b x)} \]

[In]

Int[(x^2*Sin[c + d*x])/(a + b*x)^3,x]

[Out]

-1/2*(a^2*d*Cos[c + d*x])/(b^4*(a + b*x)) - (2*a*d*Cos[c - (a*d)/b]*CosIntegral[(a*d)/b + d*x])/b^4 + (CosInte
gral[(a*d)/b + d*x]*Sin[c - (a*d)/b])/b^3 - (a^2*d^2*CosIntegral[(a*d)/b + d*x]*Sin[c - (a*d)/b])/(2*b^5) - (a
^2*Sin[c + d*x])/(2*b^3*(a + b*x)^2) + (2*a*Sin[c + d*x])/(b^3*(a + b*x)) + (Cos[c - (a*d)/b]*SinIntegral[(a*d
)/b + d*x])/b^3 - (a^2*d^2*Cos[c - (a*d)/b]*SinIntegral[(a*d)/b + d*x])/(2*b^5) + (2*a*d*Sin[c - (a*d)/b]*SinI
ntegral[(a*d)/b + d*x])/b^4

Rule 3378

Int[((c_.) + (d_.)*(x_))^(m_)*sin[(e_.) + (f_.)*(x_)], x_Symbol] :> Simp[(c + d*x)^(m + 1)*(Sin[e + f*x]/(d*(m
+ 1))), x] - Dist[f/(d*(m + 1)), Int[(c + d*x)^(m + 1)*Cos[e + f*x], x], x] /; FreeQ[{c, d, e, f}, x] && LtQ[
m, -1]

Rule 3380

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[SinIntegral[e + f*x]/d, x] /; FreeQ[{c, d,
e, f}, x] && EqQ[d*e - c*f, 0]

Rule 3383

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[CosIntegral[e - Pi/2 + f*x]/d, x] /; FreeQ
[{c, d, e, f}, x] && EqQ[d*(e - Pi/2) - c*f, 0]

Rule 3384

Int[sin[(e_.) + (f_.)*(x_)]/((c_.) + (d_.)*(x_)), x_Symbol] :> Dist[Cos[(d*e - c*f)/d], Int[Sin[c*(f/d) + f*x]
/(c + d*x), x], x] + Dist[Sin[(d*e - c*f)/d], Int[Cos[c*(f/d) + f*x]/(c + d*x), x], x] /; FreeQ[{c, d, e, f},
x] && NeQ[d*e - c*f, 0]

Rule 6874

Int[u_, x_Symbol] :> With[{v = ExpandIntegrand[u, x]}, Int[v, x] /; SumQ[v]]

Rubi steps \begin{align*} \text {integral}& = \int \left (\frac {a^2 \sin (c+d x)}{b^2 (a+b x)^3}-\frac {2 a \sin (c+d x)}{b^2 (a+b x)^2}+\frac {\sin (c+d x)}{b^2 (a+b x)}\right ) \, dx \\ & = \frac {\int \frac {\sin (c+d x)}{a+b x} \, dx}{b^2}-\frac {(2 a) \int \frac {\sin (c+d x)}{(a+b x)^2} \, dx}{b^2}+\frac {a^2 \int \frac {\sin (c+d x)}{(a+b x)^3} \, dx}{b^2} \\ & = -\frac {a^2 \sin (c+d x)}{2 b^3 (a+b x)^2}+\frac {2 a \sin (c+d x)}{b^3 (a+b x)}-\frac {(2 a d) \int \frac {\cos (c+d x)}{a+b x} \, dx}{b^3}+\frac {\left (a^2 d\right ) \int \frac {\cos (c+d x)}{(a+b x)^2} \, dx}{2 b^3}+\frac {\cos \left (c-\frac {a d}{b}\right ) \int \frac {\sin \left (\frac {a d}{b}+d x\right )}{a+b x} \, dx}{b^2}+\frac {\sin \left (c-\frac {a d}{b}\right ) \int \frac {\cos \left (\frac {a d}{b}+d x\right )}{a+b x} \, dx}{b^2} \\ & = -\frac {a^2 d \cos (c+d x)}{2 b^4 (a+b x)}+\frac {\operatorname {CosIntegral}\left (\frac {a d}{b}+d x\right ) \sin \left (c-\frac {a d}{b}\right )}{b^3}-\frac {a^2 \sin (c+d x)}{2 b^3 (a+b x)^2}+\frac {2 a \sin (c+d x)}{b^3 (a+b x)}+\frac {\cos \left (c-\frac {a d}{b}\right ) \text {Si}\left (\frac {a d}{b}+d x\right )}{b^3}-\frac {\left (a^2 d^2\right ) \int \frac {\sin (c+d x)}{a+b x} \, dx}{2 b^4}-\frac {\left (2 a d \cos \left (c-\frac {a d}{b}\right )\right ) \int \frac {\cos \left (\frac {a d}{b}+d x\right )}{a+b x} \, dx}{b^3}+\frac {\left (2 a d \sin \left (c-\frac {a d}{b}\right )\right ) \int \frac {\sin \left (\frac {a d}{b}+d x\right )}{a+b x} \, dx}{b^3} \\ & = -\frac {a^2 d \cos (c+d x)}{2 b^4 (a+b x)}-\frac {2 a d \cos \left (c-\frac {a d}{b}\right ) \operatorname {CosIntegral}\left (\frac {a d}{b}+d x\right )}{b^4}+\frac {\operatorname {CosIntegral}\left (\frac {a d}{b}+d x\right ) \sin \left (c-\frac {a d}{b}\right )}{b^3}-\frac {a^2 \sin (c+d x)}{2 b^3 (a+b x)^2}+\frac {2 a \sin (c+d x)}{b^3 (a+b x)}+\frac {\cos \left (c-\frac {a d}{b}\right ) \text {Si}\left (\frac {a d}{b}+d x\right )}{b^3}+\frac {2 a d \sin \left (c-\frac {a d}{b}\right ) \text {Si}\left (\frac {a d}{b}+d x\right )}{b^4}-\frac {\left (a^2 d^2 \cos \left (c-\frac {a d}{b}\right )\right ) \int \frac {\sin \left (\frac {a d}{b}+d x\right )}{a+b x} \, dx}{2 b^4}-\frac {\left (a^2 d^2 \sin \left (c-\frac {a d}{b}\right )\right ) \int \frac {\cos \left (\frac {a d}{b}+d x\right )}{a+b x} \, dx}{2 b^4} \\ & = -\frac {a^2 d \cos (c+d x)}{2 b^4 (a+b x)}-\frac {2 a d \cos \left (c-\frac {a d}{b}\right ) \operatorname {CosIntegral}\left (\frac {a d}{b}+d x\right )}{b^4}+\frac {\operatorname {CosIntegral}\left (\frac {a d}{b}+d x\right ) \sin \left (c-\frac {a d}{b}\right )}{b^3}-\frac {a^2 d^2 \operatorname {CosIntegral}\left (\frac {a d}{b}+d x\right ) \sin \left (c-\frac {a d}{b}\right )}{2 b^5}-\frac {a^2 \sin (c+d x)}{2 b^3 (a+b x)^2}+\frac {2 a \sin (c+d x)}{b^3 (a+b x)}+\frac {\cos \left (c-\frac {a d}{b}\right ) \text {Si}\left (\frac {a d}{b}+d x\right )}{b^3}-\frac {a^2 d^2 \cos \left (c-\frac {a d}{b}\right ) \text {Si}\left (\frac {a d}{b}+d x\right )}{2 b^5}+\frac {2 a d \sin \left (c-\frac {a d}{b}\right ) \text {Si}\left (\frac {a d}{b}+d x\right )}{b^4} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.85 (sec) , antiderivative size = 154, normalized size of antiderivative = 0.64 \[ \int \frac {x^2 \sin (c+d x)}{(a+b x)^3} \, dx=-\frac {-\operatorname {CosIntegral}\left (d \left (\frac {a}{b}+x\right )\right ) \left (-4 a b d \cos \left (c-\frac {a d}{b}\right )+\left (2 b^2-a^2 d^2\right ) \sin \left (c-\frac {a d}{b}\right )\right )+\frac {a b (a d (a+b x) \cos (c+d x)-b (3 a+4 b x) \sin (c+d x))}{(a+b x)^2}+\left (\left (-2 b^2+a^2 d^2\right ) \cos \left (c-\frac {a d}{b}\right )-4 a b d \sin \left (c-\frac {a d}{b}\right )\right ) \text {Si}\left (d \left (\frac {a}{b}+x\right )\right )}{2 b^5} \]

[In]

Integrate[(x^2*Sin[c + d*x])/(a + b*x)^3,x]

[Out]

-1/2*(-(CosIntegral[d*(a/b + x)]*(-4*a*b*d*Cos[c - (a*d)/b] + (2*b^2 - a^2*d^2)*Sin[c - (a*d)/b])) + (a*b*(a*d
*(a + b*x)*Cos[c + d*x] - b*(3*a + 4*b*x)*Sin[c + d*x]))/(a + b*x)^2 + ((-2*b^2 + a^2*d^2)*Cos[c - (a*d)/b] -
4*a*b*d*Sin[c - (a*d)/b])*SinIntegral[d*(a/b + x)])/b^5

Maple [C] (verified)

Result contains complex when optimal does not.

Time = 0.39 (sec) , antiderivative size = 621, normalized size of antiderivative = 2.58

method	result	size
risch	\(-\frac {i \left (2 i a^{2} b^{3} d^{4} x^{3}+6 i a^{3} b^{2} d^{4} x^{2}+6 i a^{4} b \,d^{4} x +2 i a^{5} d^{4}\right ) \cos \left (d x +c \right )}{4 b^{4} d \left (b x +a \right )^{2} \left (-d^{2} x^{2} b^{2}-2 a b \,d^{2} x -d^{2} a^{2}\right )}-\frac {\left (8 a \,b^{3} d^{3} x^{3}+22 a^{2} b^{2} d^{3} x^{2}+20 a^{3} b \,d^{3} x +6 a^{4} d^{3}\right ) \sin \left (d x +c \right )}{4 b^{3} d \left (b x +a \right )^{2} \left (-d^{2} x^{2} b^{2}-2 a b \,d^{2} x -d^{2} a^{2}\right )}+\frac {i \cos \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (\frac {i d \left (b x +a \right )}{b}\right ) a^{2} d^{2}}{4 b^{5}}-\frac {i \cos \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (-\frac {i d \left (b x +a \right )}{b}\right ) a^{2} d^{2}}{4 b^{5}}-\frac {i \cos \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (\frac {i d \left (b x +a \right )}{b}\right )}{2 b^{3}}+\frac {\cos \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (\frac {i d \left (b x +a \right )}{b}\right ) a d}{b^{4}}+\frac {i \cos \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (-\frac {i d \left (b x +a \right )}{b}\right )}{2 b^{3}}+\frac {\cos \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (-\frac {i d \left (b x +a \right )}{b}\right ) a d}{b^{4}}-\frac {\sin \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (\frac {i d \left (b x +a \right )}{b}\right ) a^{2} d^{2}}{4 b^{5}}-\frac {\sin \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (-\frac {i d \left (b x +a \right )}{b}\right ) a^{2} d^{2}}{4 b^{5}}+\frac {\sin \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (\frac {i d \left (b x +a \right )}{b}\right )}{2 b^{3}}+\frac {i \sin \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (\frac {i d \left (b x +a \right )}{b}\right ) a d}{b^{4}}+\frac {\sin \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (-\frac {i d \left (b x +a \right )}{b}\right )}{2 b^{3}}-\frac {i \sin \left (\frac {d a -c b}{b}\right ) \operatorname {Ei}_{1}\left (-\frac {i d \left (b x +a \right )}{b}\right ) a d}{b^{4}}\)	\(621\)
derivativedivides	\(\frac {d^{3} c^{2} \left (-\frac {\sin \left (d x +c \right )}{2 \left (d a -c b +b \left (d x +c \right )\right )^{2} b}+\frac {-\frac {\cos \left (d x +c \right )}{\left (d a -c b +b \left (d x +c \right )\right ) b}-\frac {\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}-\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}}{b}}{2 b}\right )+\frac {2 d^{3} \left (d a -c b \right ) c \left (-\frac {\sin \left (d x +c \right )}{2 \left (d a -c b +b \left (d x +c \right )\right )^{2} b}+\frac {-\frac {\cos \left (d x +c \right )}{\left (d a -c b +b \left (d x +c \right )\right ) b}-\frac {\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}-\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}}{b}}{2 b}\right )}{b}-\frac {2 d^{3} c \left (-\frac {\sin \left (d x +c \right )}{\left (d a -c b +b \left (d x +c \right )\right ) b}+\frac {\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}+\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}}{b}\right )}{b}+\frac {d^{3} \left (\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}-\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}\right )}{b^{2}}+\frac {d^{3} \left (d a -c b \right )^{2} \left (-\frac {\sin \left (d x +c \right )}{2 \left (d a -c b +b \left (d x +c \right )\right )^{2} b}+\frac {-\frac {\cos \left (d x +c \right )}{\left (d a -c b +b \left (d x +c \right )\right ) b}-\frac {\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}-\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}}{b}}{2 b}\right )}{b^{2}}-\frac {2 d^{3} \left (d a -c b \right ) \left (-\frac {\sin \left (d x +c \right )}{\left (d a -c b +b \left (d x +c \right )\right ) b}+\frac {\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}+\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}}{b}\right )}{b^{2}}}{d^{3}}\)	\(779\)
default	\(\frac {d^{3} c^{2} \left (-\frac {\sin \left (d x +c \right )}{2 \left (d a -c b +b \left (d x +c \right )\right )^{2} b}+\frac {-\frac {\cos \left (d x +c \right )}{\left (d a -c b +b \left (d x +c \right )\right ) b}-\frac {\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}-\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}}{b}}{2 b}\right )+\frac {2 d^{3} \left (d a -c b \right ) c \left (-\frac {\sin \left (d x +c \right )}{2 \left (d a -c b +b \left (d x +c \right )\right )^{2} b}+\frac {-\frac {\cos \left (d x +c \right )}{\left (d a -c b +b \left (d x +c \right )\right ) b}-\frac {\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}-\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}}{b}}{2 b}\right )}{b}-\frac {2 d^{3} c \left (-\frac {\sin \left (d x +c \right )}{\left (d a -c b +b \left (d x +c \right )\right ) b}+\frac {\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}+\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}}{b}\right )}{b}+\frac {d^{3} \left (\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}-\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}\right )}{b^{2}}+\frac {d^{3} \left (d a -c b \right )^{2} \left (-\frac {\sin \left (d x +c \right )}{2 \left (d a -c b +b \left (d x +c \right )\right )^{2} b}+\frac {-\frac {\cos \left (d x +c \right )}{\left (d a -c b +b \left (d x +c \right )\right ) b}-\frac {\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}-\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}}{b}}{2 b}\right )}{b^{2}}-\frac {2 d^{3} \left (d a -c b \right ) \left (-\frac {\sin \left (d x +c \right )}{\left (d a -c b +b \left (d x +c \right )\right ) b}+\frac {\frac {\operatorname {Si}\left (d x +c +\frac {d a -c b}{b}\right ) \sin \left (\frac {d a -c b}{b}\right )}{b}+\frac {\operatorname {Ci}\left (d x +c +\frac {d a -c b}{b}\right ) \cos \left (\frac {d a -c b}{b}\right )}{b}}{b}\right )}{b^{2}}}{d^{3}}\)	\(779\)

[In]

int(x^2*sin(d*x+c)/(b*x+a)^3,x,method=_RETURNVERBOSE)

[Out]

-1/4*I/b^4/d*(2*I*a^2*b^3*d^4*x^3+6*I*a^3*b^2*d^4*x^2+6*I*a^4*b*d^4*x+2*I*a^5*d^4)/(b*x+a)^2/(-b^2*d^2*x^2-2*a
*b*d^2*x-a^2*d^2)*cos(d*x+c)-1/4/b^3/d*(8*a*b^3*d^3*x^3+22*a^2*b^2*d^3*x^2+20*a^3*b*d^3*x+6*a^4*d^3)/(b*x+a)^2
/(-b^2*d^2*x^2-2*a*b*d^2*x-a^2*d^2)*sin(d*x+c)+1/4*I/b^5*cos((a*d-b*c)/b)*Ei(1,I*d*(b*x+a)/b)*a^2*d^2-1/4*I/b^
5*cos((a*d-b*c)/b)*Ei(1,-I*d*(b*x+a)/b)*a^2*d^2-1/2*I/b^3*cos((a*d-b*c)/b)*Ei(1,I*d*(b*x+a)/b)+1/b^4*cos((a*d-
b*c)/b)*Ei(1,I*d*(b*x+a)/b)*a*d+1/2*I/b^3*cos((a*d-b*c)/b)*Ei(1,-I*d*(b*x+a)/b)+1/b^4*cos((a*d-b*c)/b)*Ei(1,-I
*d*(b*x+a)/b)*a*d-1/4/b^5*sin((a*d-b*c)/b)*Ei(1,I*d*(b*x+a)/b)*a^2*d^2-1/4/b^5*sin((a*d-b*c)/b)*Ei(1,-I*d*(b*x
+a)/b)*a^2*d^2+1/2/b^3*sin((a*d-b*c)/b)*Ei(1,I*d*(b*x+a)/b)+I/b^4*sin((a*d-b*c)/b)*Ei(1,I*d*(b*x+a)/b)*a*d+1/2
/b^3*sin((a*d-b*c)/b)*Ei(1,-I*d*(b*x+a)/b)-I/b^4*sin((a*d-b*c)/b)*Ei(1,-I*d*(b*x+a)/b)*a*d

Fricas [A] (verification not implemented)

none

Time = 0.30 (sec) , antiderivative size = 325, normalized size of antiderivative = 1.35 \[ \int \frac {x^2 \sin (c+d x)}{(a+b x)^3} \, dx=-\frac {{\left (a^{2} b^{2} d x + a^{3} b d\right )} \cos \left (d x + c\right ) + {\left (4 \, {\left (a b^{3} d x^{2} + 2 \, a^{2} b^{2} d x + a^{3} b d\right )} \operatorname {Ci}\left (\frac {b d x + a d}{b}\right ) + {\left (a^{4} d^{2} - 2 \, a^{2} b^{2} + {\left (a^{2} b^{2} d^{2} - 2 \, b^{4}\right )} x^{2} + 2 \, {\left (a^{3} b d^{2} - 2 \, a b^{3}\right )} x\right )} \operatorname {Si}\left (\frac {b d x + a d}{b}\right )\right )} \cos \left (-\frac {b c - a d}{b}\right ) - {\left (4 \, a b^{3} x + 3 \, a^{2} b^{2}\right )} \sin \left (d x + c\right ) - {\left ({\left (a^{4} d^{2} - 2 \, a^{2} b^{2} + {\left (a^{2} b^{2} d^{2} - 2 \, b^{4}\right )} x^{2} + 2 \, {\left (a^{3} b d^{2} - 2 \, a b^{3}\right )} x\right )} \operatorname {Ci}\left (\frac {b d x + a d}{b}\right ) - 4 \, {\left (a b^{3} d x^{2} + 2 \, a^{2} b^{2} d x + a^{3} b d\right )} \operatorname {Si}\left (\frac {b d x + a d}{b}\right )\right )} \sin \left (-\frac {b c - a d}{b}\right )}{2 \, {\left (b^{7} x^{2} + 2 \, a b^{6} x + a^{2} b^{5}\right )}} \]

[In]

integrate(x^2*sin(d*x+c)/(b*x+a)^3,x, algorithm="fricas")

[Out]

-1/2*((a^2*b^2*d*x + a^3*b*d)*cos(d*x + c) + (4*(a*b^3*d*x^2 + 2*a^2*b^2*d*x + a^3*b*d)*cos_integral((b*d*x +
a*d)/b) + (a^4*d^2 - 2*a^2*b^2 + (a^2*b^2*d^2 - 2*b^4)*x^2 + 2*(a^3*b*d^2 - 2*a*b^3)*x)*sin_integral((b*d*x +
a*d)/b))*cos(-(b*c - a*d)/b) - (4*a*b^3*x + 3*a^2*b^2)*sin(d*x + c) - ((a^4*d^2 - 2*a^2*b^2 + (a^2*b^2*d^2 - 2
*b^4)*x^2 + 2*(a^3*b*d^2 - 2*a*b^3)*x)*cos_integral((b*d*x + a*d)/b) - 4*(a*b^3*d*x^2 + 2*a^2*b^2*d*x + a^3*b*
d)*sin_integral((b*d*x + a*d)/b))*sin(-(b*c - a*d)/b))/(b^7*x^2 + 2*a*b^6*x + a^2*b^5)

Sympy [F]

\[ \int \frac {x^2 \sin (c+d x)}{(a+b x)^3} \, dx=\int \frac {x^{2} \sin {\left (c + d x \right )}}{\left (a + b x\right )^{3}}\, dx \]

[In]

integrate(x**2*sin(d*x+c)/(b*x+a)**3,x)

[Out]

Integral(x**2*sin(c + d*x)/(a + b*x)**3, x)

Maxima [F]

\[ \int \frac {x^2 \sin (c+d x)}{(a+b x)^3} \, dx=\int { \frac {x^{2} \sin \left (d x + c\right )}{{\left (b x + a\right )}^{3}} \,d x } \]

[In]

integrate(x^2*sin(d*x+c)/(b*x+a)^3,x, algorithm="maxima")

[Out]

-1/2*((b*cos(c)^2 + b*sin(c)^2)*d*x^2*cos(d*x + c) + ((a*(I*exp_integral_e(4, (I*b*d*x + I*a*d)/b) - I*exp_int
egral_e(4, -(I*b*d*x + I*a*d)/b))*cos(c)^2 + a*(I*exp_integral_e(4, (I*b*d*x + I*a*d)/b) - I*exp_integral_e(4,
-(I*b*d*x + I*a*d)/b))*sin(c)^2)*cos(-(b*c - a*d)/b) - (a*(exp_integral_e(4, (I*b*d*x + I*a*d)/b) + exp_integ
ral_e(4, -(I*b*d*x + I*a*d)/b))*cos(c)^2 + a*(exp_integral_e(4, (I*b*d*x + I*a*d)/b) + exp_integral_e(4, -(I*b
*d*x + I*a*d)/b))*sin(c)^2)*sin(-(b*c - a*d)/b))*cos(d*x + c)^2 + (b*cos(c)^2 + b*sin(c)^2)*x*sin(d*x + c) + (
(a*(I*exp_integral_e(4, (I*b*d*x + I*a*d)/b) - I*exp_integral_e(4, -(I*b*d*x + I*a*d)/b))*cos(c)^2 + a*(I*exp_
integral_e(4, (I*b*d*x + I*a*d)/b) - I*exp_integral_e(4, -(I*b*d*x + I*a*d)/b))*sin(c)^2)*cos(-(b*c - a*d)/b)
- (a*(exp_integral_e(4, (I*b*d*x + I*a*d)/b) + exp_integral_e(4, -(I*b*d*x + I*a*d)/b))*cos(c)^2 + a*(exp_inte
gral_e(4, (I*b*d*x + I*a*d)/b) + exp_integral_e(4, -(I*b*d*x + I*a*d)/b))*sin(c)^2)*sin(-(b*c - a*d)/b))*sin(d
*x + c)^2 + ((b*d*x^2*cos(c) - b*x*sin(c))*cos(d*x + c)^2 + (b*d*x^2*cos(c) - b*x*sin(c))*sin(d*x + c)^2)*cos(
d*x + 2*c) - 6*(((a*b^4*cos(c)^2 + a*b^4*sin(c)^2)*d^3*x^3 + 3*(a^2*b^3*cos(c)^2 + a^2*b^3*sin(c)^2)*d^3*x^2 +
3*(a^3*b^2*cos(c)^2 + a^3*b^2*sin(c)^2)*d^3*x + (a^4*b*cos(c)^2 + a^4*b*sin(c)^2)*d^3)*cos(d*x + c)^2 + ((a*b
^4*cos(c)^2 + a*b^4*sin(c)^2)*d^3*x^3 + 3*(a^2*b^3*cos(c)^2 + a^2*b^3*sin(c)^2)*d^3*x^2 + 3*(a^3*b^2*cos(c)^2
+ a^3*b^2*sin(c)^2)*d^3*x + (a^4*b*cos(c)^2 + a^4*b*sin(c)^2)*d^3)*sin(d*x + c)^2)*integrate(1/2*x*cos(d*x + c
)/(b^4*d^2*x^4 + 4*a*b^3*d^2*x^3 + 6*a^2*b^2*d^2*x^2 + 4*a^3*b*d^2*x + a^4*d^2), x) - 6*(((a*b^4*cos(c)^2 + a*
b^4*sin(c)^2)*d^3*x^3 + 3*(a^2*b^3*cos(c)^2 + a^2*b^3*sin(c)^2)*d^3*x^2 + 3*(a^3*b^2*cos(c)^2 + a^3*b^2*sin(c)
^2)*d^3*x + (a^4*b*cos(c)^2 + a^4*b*sin(c)^2)*d^3)*cos(d*x + c)^2 + ((a*b^4*cos(c)^2 + a*b^4*sin(c)^2)*d^3*x^3
+ 3*(a^2*b^3*cos(c)^2 + a^2*b^3*sin(c)^2)*d^3*x^2 + 3*(a^3*b^2*cos(c)^2 + a^3*b^2*sin(c)^2)*d^3*x + (a^4*b*co
s(c)^2 + a^4*b*sin(c)^2)*d^3)*sin(d*x + c)^2)*integrate(1/2*x*cos(d*x + c)/((b^4*d^2*x^4 + 4*a*b^3*d^2*x^3 + 6
*a^2*b^2*d^2*x^2 + 4*a^3*b*d^2*x + a^4*d^2)*cos(d*x + c)^2 + (b^4*d^2*x^4 + 4*a*b^3*d^2*x^3 + 6*a^2*b^2*d^2*x^
2 + 4*a^3*b*d^2*x + a^4*d^2)*sin(d*x + c)^2), x) + 4*(((b^5*cos(c)^2 + b^5*sin(c)^2)*d^2*x^3 + 3*(a*b^4*cos(c)
^2 + a*b^4*sin(c)^2)*d^2*x^2 + 3*(a^2*b^3*cos(c)^2 + a^2*b^3*sin(c)^2)*d^2*x + (a^3*b^2*cos(c)^2 + a^3*b^2*sin
(c)^2)*d^2)*cos(d*x + c)^2 + ((b^5*cos(c)^2 + b^5*sin(c)^2)*d^2*x^3 + 3*(a*b^4*cos(c)^2 + a*b^4*sin(c)^2)*d^2*
x^2 + 3*(a^2*b^3*cos(c)^2 + a^2*b^3*sin(c)^2)*d^2*x + (a^3*b^2*cos(c)^2 + a^3*b^2*sin(c)^2)*d^2)*sin(d*x + c)^
2)*integrate(1/2*x*sin(d*x + c)/(b^4*d^2*x^4 + 4*a*b^3*d^2*x^3 + 6*a^2*b^2*d^2*x^2 + 4*a^3*b*d^2*x + a^4*d^2),
x) + 4*(((b^5*cos(c)^2 + b^5*sin(c)^2)*d^2*x^3 + 3*(a*b^4*cos(c)^2 + a*b^4*sin(c)^2)*d^2*x^2 + 3*(a^2*b^3*cos
(c)^2 + a^2*b^3*sin(c)^2)*d^2*x + (a^3*b^2*cos(c)^2 + a^3*b^2*sin(c)^2)*d^2)*cos(d*x + c)^2 + ((b^5*cos(c)^2 +
b^5*sin(c)^2)*d^2*x^3 + 3*(a*b^4*cos(c)^2 + a*b^4*sin(c)^2)*d^2*x^2 + 3*(a^2*b^3*cos(c)^2 + a^2*b^3*sin(c)^2)
*d^2*x + (a^3*b^2*cos(c)^2 + a^3*b^2*sin(c)^2)*d^2)*sin(d*x + c)^2)*integrate(1/2*x*sin(d*x + c)/((b^4*d^2*x^4
+ 4*a*b^3*d^2*x^3 + 6*a^2*b^2*d^2*x^2 + 4*a^3*b*d^2*x + a^4*d^2)*cos(d*x + c)^2 + (b^4*d^2*x^4 + 4*a*b^3*d^2*
x^3 + 6*a^2*b^2*d^2*x^2 + 4*a^3*b*d^2*x + a^4*d^2)*sin(d*x + c)^2), x) + ((b*d*x^2*sin(c) + b*x*cos(c))*cos(d*
x + c)^2 + (b*d*x^2*sin(c) + b*x*cos(c))*sin(d*x + c)^2)*sin(d*x + 2*c))/(((b^4*cos(c)^2 + b^4*sin(c)^2)*d^2*x
^3 + 3*(a*b^3*cos(c)^2 + a*b^3*sin(c)^2)*d^2*x^2 + 3*(a^2*b^2*cos(c)^2 + a^2*b^2*sin(c)^2)*d^2*x + (a^3*b*cos(
c)^2 + a^3*b*sin(c)^2)*d^2)*cos(d*x + c)^2 + ((b^4*cos(c)^2 + b^4*sin(c)^2)*d^2*x^3 + 3*(a*b^3*cos(c)^2 + a*b^
3*sin(c)^2)*d^2*x^2 + 3*(a^2*b^2*cos(c)^2 + a^2*b^2*sin(c)^2)*d^2*x + (a^3*b*cos(c)^2 + a^3*b*sin(c)^2)*d^2)*s
in(d*x + c)^2)

Giac [C] (verification not implemented)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 0.62 (sec) , antiderivative size = 15410, normalized size of antiderivative = 63.94 \[ \int \frac {x^2 \sin (c+d x)}{(a+b x)^3} \, dx=\text {Too large to display} \]

[In]

integrate(x^2*sin(d*x+c)/(b*x+a)^3,x, algorithm="giac")

[Out]

-1/4*(a^2*b^2*d^2*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - a^2*
b^2*d^2*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*a^2*b^2*d^2
*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*a^2*b^2*d^2*x^2*real_part(
cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) + 2*a^2*b^2*d^2*x^2*real_part(cos_integr
al(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) - 2*a^2*b^2*d^2*x^2*real_part(cos_integral(d*x +
a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 - 2*a^2*b^2*d^2*x^2*real_part(cos_integral(-d*x - a*d/b))*t
an(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 + 2*a^3*b*d^2*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*
tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 2*a^3*b*d^2*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^
2*tan(1/2*a*d/b)^2 + 4*a*b^3*d*x^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*
d/b)^2 + 4*a*b^3*d*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 4*
a^3*b*d^2*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - a^2*b^2*d^2*x^2*imag_
part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 + a^2*b^2*d^2*x^2*imag_part(cos_integral(-d*x - a*
d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*a^2*b^2*d^2*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^
2 + 4*a^2*b^2*d^2*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) - 4*a^2*b^
2*d^2*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) + 8*a^2*b^2*d^2*x^2*s
in_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) - 8*a*b^3*d*x^2*imag_part(cos_integral(d
*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) + 8*a*b^3*d*x^2*imag_part(cos_integral(-d*x - a*d/b))*
tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) + 4*a^3*b*d^2*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2
*tan(1/2*c)^2*tan(1/2*a*d/b) + 4*a^3*b*d^2*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2
*tan(1/2*a*d/b) - 16*a*b^3*d*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) - a^
2*b^2*d^2*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + a^2*b^2*d^2*x^2*imag_part
(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 - 2*a^2*b^2*d^2*x^2*sin_integral((b*d*x + a*d)/b)
*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + 8*a*b^3*d*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c
)*tan(1/2*a*d/b)^2 - 8*a*b^3*d*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d
/b)^2 - 4*a^3*b*d^2*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 - 4*a^3*
b*d^2*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 + 16*a*b^3*d*x^2*sin_
integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 + a^2*b^2*d^2*x^2*imag_part(cos_integral(
d*x + a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - a^2*b^2*d^2*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)
^2*tan(1/2*a*d/b)^2 + 2*a^2*b^2*d^2*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + a^4*d^2*
imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 2*b^4*x^2*imag_part(cos_in
tegral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - a^4*d^2*imag_part(cos_integral(-d*x - a*d/
b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*b^4*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x
)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 8*a^2*b^2*d*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*
c)^2*tan(1/2*a*d/b)^2 + 8*a^2*b^2*d*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/
2*a*d/b)^2 + 2*a^4*d^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 4*b^4*x^2*
sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*a^2*b^2*d^2*x^2*real_part(cos_i
ntegral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c) + 2*a^2*b^2*d^2*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(
1/2*d*x)^2*tan(1/2*c) - 2*a^3*b*d^2*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*a^3
*b*d^2*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 4*a*b^3*d*x^2*real_part(cos_integ
ral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 4*a*b^3*d*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*
x)^2*tan(1/2*c)^2 - 4*a^3*b*d^2*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*a^2*b^2*d^2*x^
2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) - 2*a^2*b^2*d^2*x^2*real_part(cos_integra
l(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) + 8*a^3*b*d^2*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*
x)^2*tan(1/2*c)*tan(1/2*a*d/b) - 8*a^3*b*d^2*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)
*tan(1/2*a*d/b) + 16*a*b^3*d*x^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)
+ 16*a*b^3*d*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) + 16*a^3*b*d^
2*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) + 2*a^2*b^2*d^2*x^2*real_part(cos_i
ntegral(d*x + a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) + 2*a^2*b^2*d^2*x^2*real_part(cos_integral(-d*x - a*d/b))*ta
n(1/2*c)^2*tan(1/2*a*d/b) - 16*a^2*b^2*d*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*ta
n(1/2*a*d/b) + 16*a^2*b^2*d*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)
+ 2*a^4*d^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) - 4*b^4*x^2*real_
part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) + 2*a^4*d^2*real_part(cos_integral(
-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) - 4*b^4*x^2*real_part(cos_integral(-d*x - a*d/b))*ta
n(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) - 32*a^2*b^2*d*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2
*c)^2*tan(1/2*a*d/b) - 2*a^3*b*d^2*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + 2*
a^3*b*d^2*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 - 4*a*b^3*d*x^2*real_part(co
s_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 - 4*a*b^3*d*x^2*real_part(cos_integral(-d*x - a*d/b))
*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 - 4*a^3*b*d^2*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2
- 2*a^2*b^2*d^2*x^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 - 2*a^2*b^2*d^2*x^2*real
_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 + 16*a^2*b^2*d*x*imag_part(cos_integral(d*x + a*
d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 - 16*a^2*b^2*d*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/
2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 - 2*a^4*d^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c
)*tan(1/2*a*d/b)^2 + 4*b^4*x^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2
- 2*a^4*d^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 + 4*b^4*x^2*real
_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 + 32*a^2*b^2*d*x*sin_integral((b*
d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 + 2*a^3*b*d^2*x*imag_part(cos_integral(d*x + a*d/b))*
tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 2*a^3*b*d^2*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b
)^2 + 4*a*b^3*d*x^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 4*a*b^3*d*x^2*real_pa
rt(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 4*a^3*b*d^2*x*sin_integral((b*d*x + a*d)/b)*tan
(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*a^2*b^2*d*x*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 4*a*b^3*x*imag_part(
cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 4*a*b^3*x*imag_part(cos_integral(-d*
x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 4*a^3*b*d*real_part(cos_integral(d*x + a*d/b))*tan(
1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 4*a^3*b*d*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(
1/2*c)^2*tan(1/2*a*d/b)^2 - 8*a*b^3*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)
^2 + a^2*b^2*d^2*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2 - a^2*b^2*d^2*x^2*imag_part(cos_integ
ral(-d*x - a*d/b))*tan(1/2*d*x)^2 + 2*a^2*b^2*d^2*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2 - 8*a*b^3*d
*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c) + 8*a*b^3*d*x^2*imag_part(cos_integral(-d*
x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c) + 4*a^3*b*d^2*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(
1/2*c) + 4*a^3*b*d^2*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c) - 16*a*b^3*d*x^2*sin_in
tegral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c) - a^2*b^2*d^2*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1
/2*c)^2 + a^2*b^2*d^2*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2 - 2*a^2*b^2*d^2*x^2*sin_integral(
(b*d*x + a*d)/b)*tan(1/2*c)^2 - a^4*d^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*b
^4*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 + a^4*d^2*imag_part(cos_integral(-d*x
- a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*b^4*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2
*c)^2 - 8*a^2*b^2*d*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 8*a^2*b^2*d*x*real_pa
rt(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*a^4*d^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d
*x)^2*tan(1/2*c)^2 + 4*b^4*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^2 + 8*a*b^3*d*x^2*imag_
part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) - 8*a*b^3*d*x^2*imag_part(cos_integral(-d*x - a*
d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) - 4*a^3*b*d^2*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/
2*a*d/b) - 4*a^3*b*d^2*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) + 16*a*b^3*d*x^2*
sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*a*d/b) + 4*a^2*b^2*d^2*x^2*imag_part(cos_integral(d*x + a
*d/b))*tan(1/2*c)*tan(1/2*a*d/b) - 4*a^2*b^2*d^2*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/2*
a*d/b) + 8*a^2*b^2*d^2*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)*tan(1/2*a*d/b) + 4*a^4*d^2*imag_part(cos_i
ntegral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) - 8*b^4*x^2*imag_part(cos_integral(d*x + a*d/b)
)*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) - 4*a^4*d^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*ta
n(1/2*c)*tan(1/2*a*d/b) + 8*b^4*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*
d/b) + 32*a^2*b^2*d*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) + 32*a^2*b
^2*d*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) + 8*a^4*d^2*sin_integral
((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) - 16*b^4*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*
d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) - 8*a*b^3*d*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2*tan(1/2*a*d
/b) + 8*a*b^3*d*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) + 4*a^3*b*d^2*x*real_par
t(cos_integral(d*x + a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) + 4*a^3*b*d^2*x*real_part(cos_integral(-d*x - a*d/b))
*tan(1/2*c)^2*tan(1/2*a*d/b) - 16*a*b^3*d*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)^2*tan(1/2*a*d/b) - 8*a^
3*b*d*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) + 8*a^3*b*d*imag_part(co
s_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) - 8*a*b^3*x*real_part(cos_integral(d*x +
a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) - 8*a*b^3*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d
*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b) - 16*a^3*b*d*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1
/2*a*d/b) - a^2*b^2*d^2*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b)^2 + a^2*b^2*d^2*x^2*imag_part(
cos_integral(-d*x - a*d/b))*tan(1/2*a*d/b)^2 - 2*a^2*b^2*d^2*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*a*d/b)^
2 - a^4*d^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + 2*b^4*x^2*imag_part(cos_int
egral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + a^4*d^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*
x)^2*tan(1/2*a*d/b)^2 - 2*b^4*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 - 8*a^
2*b^2*d*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 - 8*a^2*b^2*d*x*real_part(cos_i
ntegral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 - 2*a^4*d^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^
2*tan(1/2*a*d/b)^2 + 4*b^4*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + 8*a*b^3*d*x^2*i
mag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 - 8*a*b^3*d*x^2*imag_part(cos_integral(-d*x -
a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 - 4*a^3*b*d^2*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*
d/b)^2 - 4*a^3*b*d^2*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 + 16*a*b^3*d*x^2*sin_
integral((b*d*x + a*d)/b)*tan(1/2*c)*tan(1/2*a*d/b)^2 + 8*a^3*b*d*imag_part(cos_integral(d*x + a*d/b))*tan(1/2
*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 - 8*a^3*b*d*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c
)*tan(1/2*a*d/b)^2 + 8*a*b^3*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2
+ 8*a*b^3*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 + 16*a^3*b*d*sin
_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 + a^4*d^2*imag_part(cos_integral(d*x + a
*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 2*b^4*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2*tan(1/2*a*d
/b)^2 - a^4*d^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*b^4*x^2*imag_part(cos_
integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 8*a^2*b^2*d*x*real_part(cos_integral(d*x + a*d/b))*tan
(1/2*c)^2*tan(1/2*a*d/b)^2 + 8*a^2*b^2*d*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b)^2
+ 2*a^4*d^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 4*b^4*x^2*sin_integral((b*d*x + a*d
)/b)*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*a^3*b*d*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 2*a^2*b^2*imag_p
art(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*a^2*b^2*imag_part(cos_integral
(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 4*a^2*b^2*sin_integral((b*d*x + a*d)/b)*tan(1/2
*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*a^3*b*d^2*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2 - 2*
a^3*b*d^2*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2 + 4*a*b^3*d*x^2*real_part(cos_integral(d*x +
a*d/b))*tan(1/2*d*x)^2 + 4*a*b^3*d*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2 + 4*a^3*b*d^2*x*si
n_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2 + 2*a^2*b^2*d^2*x^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)
+ 2*a^2*b^2*d^2*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c) - 16*a^2*b^2*d*x*imag_part(cos_integral(
d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c) + 16*a^2*b^2*d*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*
tan(1/2*c) + 2*a^4*d^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c) - 4*b^4*x^2*real_part(co
s_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c) + 2*a^4*d^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d
*x)^2*tan(1/2*c) - 4*b^4*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c) - 32*a^2*b^2*d*x*
sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c) - 2*a^3*b*d^2*x*imag_part(cos_integral(d*x + a*d/b))*t
an(1/2*c)^2 + 2*a^3*b*d^2*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2 - 4*a*b^3*d*x^2*real_part(cos_i
ntegral(d*x + a*d/b))*tan(1/2*c)^2 - 4*a*b^3*d*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2 - 4*a^3*
b*d^2*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)^2 + 2*a^2*b^2*d*x*tan(1/2*d*x)^2*tan(1/2*c)^2 + 4*a*b^3*x*ima
g_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 4*a*b^3*x*imag_part(cos_integral(-d*x - a*d/b)
)*tan(1/2*d*x)^2*tan(1/2*c)^2 - 4*a^3*b*d*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 4
*a^3*b*d*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 + 8*a*b^3*x*sin_integral((b*d*x + a
*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*a^2*b^2*d^2*x^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b) - 2
*a^2*b^2*d^2*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*a*d/b) + 16*a^2*b^2*d*x*imag_part(cos_integral(
d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) - 16*a^2*b^2*d*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x
)^2*tan(1/2*a*d/b) - 2*a^4*d^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) + 4*b^4*x^2*
real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) - 2*a^4*d^2*real_part(cos_integral(-d*x - a
*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) + 4*b^4*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*
a*d/b) + 32*a^2*b^2*d*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*a*d/b) + 8*a^3*b*d^2*x*imag_part(
cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) - 8*a^3*b*d^2*x*imag_part(cos_integral(-d*x - a*d/b))*tan
(1/2*c)*tan(1/2*a*d/b) + 16*a*b^3*d*x^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) + 16*a*
b^3*d*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) + 16*a^3*b*d^2*x*sin_integral((b*d*x
+ a*d)/b)*tan(1/2*c)*tan(1/2*a*d/b) - 16*a*b^3*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*
c)*tan(1/2*a*d/b) + 16*a*b^3*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)
+ 16*a^3*b*d*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) + 16*a^3*b*d*real_p
art(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) - 32*a*b^3*x*sin_integral((b*d*x + a*
d)/b)*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) - 16*a^2*b^2*d*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c
)^2*tan(1/2*a*d/b) + 16*a^2*b^2*d*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) + 2*a^4*
d^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) - 4*b^4*x^2*real_part(cos_integral(d*x +
a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) + 2*a^4*d^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d
/b) - 4*b^4*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) - 32*a^2*b^2*d*x*sin_integra
l((b*d*x + a*d)/b)*tan(1/2*c)^2*tan(1/2*a*d/b) - 4*a^2*b^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2
*tan(1/2*c)^2*tan(1/2*a*d/b) - 4*a^2*b^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2*tan
(1/2*a*d/b) - 2*a^3*b*d^2*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b)^2 + 2*a^3*b*d^2*x*imag_part(co
s_integral(-d*x - a*d/b))*tan(1/2*a*d/b)^2 - 4*a*b^3*d*x^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b)
^2 - 4*a*b^3*d*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*a*d/b)^2 - 4*a^3*b*d^2*x*sin_integral((b*d*x
+ a*d)/b)*tan(1/2*a*d/b)^2 - 2*a^2*b^2*d*x*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + 4*a*b^3*x*imag_part(cos_integral(
d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 - 4*a*b^3*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2
*tan(1/2*a*d/b)^2 - 4*a^3*b*d*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 - 4*a^3*b*d
*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + 8*a*b^3*x*sin_integral((b*d*x + a*d)/
b)*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + 16*a^2*b^2*d*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*
d/b)^2 - 16*a^2*b^2*d*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 - 2*a^4*d^2*real_par
t(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 + 4*b^4*x^2*real_part(cos_integral(d*x + a*d/b))*tan(
1/2*c)*tan(1/2*a*d/b)^2 - 2*a^4*d^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 + 4*b^4*
x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 + 32*a^2*b^2*d*x*sin_integral((b*d*x + a
*d)/b)*tan(1/2*c)*tan(1/2*a*d/b)^2 - 8*a^2*b^2*d*x*tan(1/2*d*x)*tan(1/2*c)*tan(1/2*a*d/b)^2 + 16*a*b^3*x*tan(1
/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 + 4*a^2*b^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*
c)*tan(1/2*a*d/b)^2 + 4*a^2*b^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)
^2 - 2*a^2*b^2*d*x*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 4*a*b^3*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2
*tan(1/2*a*d/b)^2 + 4*a*b^3*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 4*a^3*b*d*
real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 4*a^3*b*d*real_part(cos_integral(-d*x - a
*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 8*a*b^3*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)^2*tan(1/2*a*d/b)^2 +
16*a*b^3*x*tan(1/2*d*x)*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + a^2*b^2*d^2*x^2*imag_part(cos_integral(d*x + a*d/b))
- a^2*b^2*d^2*x^2*imag_part(cos_integral(-d*x - a*d/b)) + 2*a^2*b^2*d^2*x^2*sin_integral((b*d*x + a*d)/b) + a^
4*d^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2 - 2*b^4*x^2*imag_part(cos_integral(d*x + a*d/b))*tan
(1/2*d*x)^2 - a^4*d^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2 + 2*b^4*x^2*imag_part(cos_integral(
-d*x - a*d/b))*tan(1/2*d*x)^2 + 8*a^2*b^2*d*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2 + 8*a^2*b^2*
d*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2 + 2*a^4*d^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x
)^2 - 4*b^4*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2 - 8*a*b^3*d*x^2*imag_part(cos_integral(d*x + a*d/
b))*tan(1/2*c) + 8*a*b^3*d*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c) + 4*a^3*b*d^2*x*real_part(cos_
integral(d*x + a*d/b))*tan(1/2*c) + 4*a^3*b*d^2*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c) - 16*a*b^3*
d*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*c) - 8*a^3*b*d*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2
*tan(1/2*c) + 8*a^3*b*d*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c) - 8*a*b^3*x*real_part(
cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c) - 8*a*b^3*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2
*d*x)^2*tan(1/2*c) - 16*a^3*b*d*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c) - a^4*d^2*imag_part(co
s_integral(d*x + a*d/b))*tan(1/2*c)^2 + 2*b^4*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2 + a^4*d^2*
imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2 - 2*b^4*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c
)^2 - 8*a^2*b^2*d*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2 - 8*a^2*b^2*d*x*real_part(cos_integral(-
d*x - a*d/b))*tan(1/2*c)^2 - 2*a^4*d^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)^2 + 4*b^4*x^2*sin_integral((b*
d*x + a*d)/b)*tan(1/2*c)^2 + 2*a^3*b*d*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*a^2*b^2*imag_part(cos_integral(d*x + a*
d/b))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*a^2*b^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)^
2 + 4*a^2*b^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)^2 + 8*a*b^3*d*x^2*imag_part(cos_integral
(d*x + a*d/b))*tan(1/2*a*d/b) - 8*a*b^3*d*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*a*d/b) - 4*a^3*b*d
^2*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b) - 4*a^3*b*d^2*x*real_part(cos_integral(-d*x - a*d/b))
*tan(1/2*a*d/b) + 16*a*b^3*d*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*a*d/b) + 8*a^3*b*d*imag_part(cos_integr
al(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) - 8*a^3*b*d*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^
2*tan(1/2*a*d/b) + 8*a*b^3*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) + 8*a*b^3*x*re
al_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) + 16*a^3*b*d*sin_integral((b*d*x + a*d)/b)*t
an(1/2*d*x)^2*tan(1/2*a*d/b) + 4*a^4*d^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) - 8*b^
4*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) - 4*a^4*d^2*imag_part(cos_integral(-d*x -
a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) + 8*b^4*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)
+ 32*a^2*b^2*d*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) + 32*a^2*b^2*d*x*real_part(co
s_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) + 8*a^4*d^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)*tan(1
/2*a*d/b) - 16*b^4*x^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)*tan(1/2*a*d/b) - 8*a^2*b^2*imag_part(cos_integ
ral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) + 8*a^2*b^2*imag_part(cos_integral(-d*x - a*d/b))*t
an(1/2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b) - 16*a^2*b^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*c)*t
an(1/2*a*d/b) - 8*a^3*b*d*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) + 8*a^3*b*d*imag_pa
rt(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) - 8*a*b^3*x*real_part(cos_integral(d*x + a*d/b))*ta
n(1/2*c)^2*tan(1/2*a*d/b) - 8*a*b^3*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) - 16*a
^3*b*d*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)^2*tan(1/2*a*d/b) - a^4*d^2*imag_part(cos_integral(d*x + a*d/b)
)*tan(1/2*a*d/b)^2 + 2*b^4*x^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b)^2 + a^4*d^2*imag_part(cos_i
ntegral(-d*x - a*d/b))*tan(1/2*a*d/b)^2 - 2*b^4*x^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*a*d/b)^2 - 8
*a^2*b^2*d*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b)^2 - 8*a^2*b^2*d*x*real_part(cos_integral(-d*x
- a*d/b))*tan(1/2*a*d/b)^2 - 2*a^4*d^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*a*d/b)^2 + 4*b^4*x^2*sin_integra
l((b*d*x + a*d)/b)*tan(1/2*a*d/b)^2 - 2*a^3*b*d*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + 2*a^2*b^2*imag_part(cos_inte
gral(d*x + a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 - 2*a^2*b^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d
*x)^2*tan(1/2*a*d/b)^2 + 4*a^2*b^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + 8*a^3*b*d*i
mag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 - 8*a^3*b*d*imag_part(cos_integral(-d*x - a*d/
b))*tan(1/2*c)*tan(1/2*a*d/b)^2 + 8*a*b^3*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 +
8*a*b^3*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 + 16*a^3*b*d*sin_integral((b*d*x
+ a*d)/b)*tan(1/2*c)*tan(1/2*a*d/b)^2 - 8*a^3*b*d*tan(1/2*d*x)*tan(1/2*c)*tan(1/2*a*d/b)^2 + 12*a^2*b^2*tan(1/
2*d*x)^2*tan(1/2*c)*tan(1/2*a*d/b)^2 - 2*a^3*b*d*tan(1/2*c)^2*tan(1/2*a*d/b)^2 - 2*a^2*b^2*imag_part(cos_integ
ral(d*x + a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*a^2*b^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2
*tan(1/2*a*d/b)^2 - 4*a^2*b^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 12*a^2*b^2*tan(1/2
*d*x)*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*a^3*b*d^2*x*imag_part(cos_integral(d*x + a*d/b)) - 2*a^3*b*d^2*x*imag_
part(cos_integral(-d*x - a*d/b)) + 4*a*b^3*d*x^2*real_part(cos_integral(d*x + a*d/b)) + 4*a*b^3*d*x^2*real_par
t(cos_integral(-d*x - a*d/b)) + 4*a^3*b*d^2*x*sin_integral((b*d*x + a*d)/b) - 2*a^2*b^2*d*x*tan(1/2*d*x)^2 - 4
*a*b^3*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2 + 4*a*b^3*x*imag_part(cos_integral(-d*x - a*d/b))
*tan(1/2*d*x)^2 + 4*a^3*b*d*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2 + 4*a^3*b*d*real_part(cos_inte
gral(-d*x - a*d/b))*tan(1/2*d*x)^2 - 8*a*b^3*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2 - 16*a^2*b^2*d*x*i
mag_part(cos_integral(d*x + a*d/b))*tan(1/2*c) + 16*a^2*b^2*d*x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*
c) + 2*a^4*d^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c) - 4*b^4*x^2*real_part(cos_integral(d*x + a*d/b)
)*tan(1/2*c) + 2*a^4*d^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c) - 4*b^4*x^2*real_part(cos_integral(-
d*x - a*d/b))*tan(1/2*c) - 32*a^2*b^2*d*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*c) - 8*a^2*b^2*d*x*tan(1/2*d*x
)*tan(1/2*c) + 16*a*b^3*x*tan(1/2*d*x)^2*tan(1/2*c) - 4*a^2*b^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d
*x)^2*tan(1/2*c) - 4*a^2*b^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*c) - 2*a^2*b^2*d*x*t
an(1/2*c)^2 + 4*a*b^3*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2 - 4*a*b^3*x*imag_part(cos_integral(-
d*x - a*d/b))*tan(1/2*c)^2 - 4*a^3*b*d*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2 - 4*a^3*b*d*real_part
(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2 + 8*a*b^3*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)^2 + 16*a*b^3*x*
tan(1/2*d*x)*tan(1/2*c)^2 + 16*a^2*b^2*d*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b) - 16*a^2*b^2*d*
x*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*a*d/b) - 2*a^4*d^2*real_part(cos_integral(d*x + a*d/b))*tan(1/
2*a*d/b) + 4*b^4*x^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b) - 2*a^4*d^2*real_part(cos_integral(-d
*x - a*d/b))*tan(1/2*a*d/b) + 4*b^4*x^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*a*d/b) + 32*a^2*b^2*d*x*
sin_integral((b*d*x + a*d)/b)*tan(1/2*a*d/b) + 4*a^2*b^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2*t
an(1/2*a*d/b) + 4*a^2*b^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^2*tan(1/2*a*d/b) - 16*a*b^3*x*ima
g_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) + 16*a*b^3*x*imag_part(cos_integral(-d*x - a*d/b))
*tan(1/2*c)*tan(1/2*a*d/b) + 16*a^3*b*d*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) + 16*a^
3*b*d*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) - 32*a*b^3*x*sin_integral((b*d*x + a*d)/
b)*tan(1/2*c)*tan(1/2*a*d/b) - 4*a^2*b^2*real_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) - 4*
a^2*b^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)^2*tan(1/2*a*d/b) + 2*a^2*b^2*d*x*tan(1/2*a*d/b)^2 + 4
*a*b^3*x*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b)^2 - 4*a*b^3*x*imag_part(cos_integral(-d*x - a*d/b
))*tan(1/2*a*d/b)^2 - 4*a^3*b*d*real_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b)^2 - 4*a^3*b*d*real_part(co
s_integral(-d*x - a*d/b))*tan(1/2*a*d/b)^2 + 8*a*b^3*x*sin_integral((b*d*x + a*d)/b)*tan(1/2*a*d/b)^2 - 16*a*b
^3*x*tan(1/2*d*x)*tan(1/2*a*d/b)^2 - 16*a*b^3*x*tan(1/2*c)*tan(1/2*a*d/b)^2 + 4*a^2*b^2*real_part(cos_integral
(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b)^2 + 4*a^2*b^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c)*tan(1/
2*a*d/b)^2 + a^4*d^2*imag_part(cos_integral(d*x + a*d/b)) - 2*b^4*x^2*imag_part(cos_integral(d*x + a*d/b)) - a
^4*d^2*imag_part(cos_integral(-d*x - a*d/b)) + 2*b^4*x^2*imag_part(cos_integral(-d*x - a*d/b)) + 8*a^2*b^2*d*x
*real_part(cos_integral(d*x + a*d/b)) + 8*a^2*b^2*d*x*real_part(cos_integral(-d*x - a*d/b)) + 2*a^4*d^2*sin_in
tegral((b*d*x + a*d)/b) - 4*b^4*x^2*sin_integral((b*d*x + a*d)/b) - 2*a^3*b*d*tan(1/2*d*x)^2 - 2*a^2*b^2*imag_
part(cos_integral(d*x + a*d/b))*tan(1/2*d*x)^2 + 2*a^2*b^2*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*d*x)^
2 - 4*a^2*b^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*d*x)^2 - 8*a^3*b*d*imag_part(cos_integral(d*x + a*d/b))*ta
n(1/2*c) + 8*a^3*b*d*imag_part(cos_integral(-d*x - a*d/b))*tan(1/2*c) - 8*a*b^3*x*real_part(cos_integral(d*x +
a*d/b))*tan(1/2*c) - 8*a*b^3*x*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c) - 16*a^3*b*d*sin_integral((b*
d*x + a*d)/b)*tan(1/2*c) - 8*a^3*b*d*tan(1/2*d*x)*tan(1/2*c) + 12*a^2*b^2*tan(1/2*d*x)^2*tan(1/2*c) - 2*a^3*b*
d*tan(1/2*c)^2 + 2*a^2*b^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)^2 - 2*a^2*b^2*imag_part(cos_integra
l(-d*x - a*d/b))*tan(1/2*c)^2 + 4*a^2*b^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)^2 + 12*a^2*b^2*tan(1/2*d*x)
*tan(1/2*c)^2 + 8*a^3*b*d*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b) - 8*a^3*b*d*imag_part(cos_integr
al(-d*x - a*d/b))*tan(1/2*a*d/b) + 8*a*b^3*x*real_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b) + 8*a*b^3*x*r
eal_part(cos_integral(-d*x - a*d/b))*tan(1/2*a*d/b) + 16*a^3*b*d*sin_integral((b*d*x + a*d)/b)*tan(1/2*a*d/b)
- 8*a^2*b^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) + 8*a^2*b^2*imag_part(cos_integral(
-d*x - a*d/b))*tan(1/2*c)*tan(1/2*a*d/b) - 16*a^2*b^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*c)*tan(1/2*a*d/b)
+ 2*a^3*b*d*tan(1/2*a*d/b)^2 + 2*a^2*b^2*imag_part(cos_integral(d*x + a*d/b))*tan(1/2*a*d/b)^2 - 2*a^2*b^2*ima
g_part(cos_integral(-d*x - a*d/b))*tan(1/2*a*d/b)^2 + 4*a^2*b^2*sin_integral((b*d*x + a*d)/b)*tan(1/2*a*d/b)^2
- 12*a^2*b^2*tan(1/2*d*x)*tan(1/2*a*d/b)^2 - 12*a^2*b^2*tan(1/2*c)*tan(1/2*a*d/b)^2 + 2*a^2*b^2*d*x - 4*a*b^3
*x*imag_part(cos_integral(d*x + a*d/b)) + 4*a*b^3*x*imag_part(cos_integral(-d*x - a*d/b)) + 4*a^3*b*d*real_par
t(cos_integral(d*x + a*d/b)) + 4*a^3*b*d*real_part(cos_integral(-d*x - a*d/b)) - 8*a*b^3*x*sin_integral((b*d*x
+ a*d)/b) - 16*a*b^3*x*tan(1/2*d*x) - 16*a*b^3*x*tan(1/2*c) - 4*a^2*b^2*real_part(cos_integral(d*x + a*d/b))*
tan(1/2*c) - 4*a^2*b^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*c) + 4*a^2*b^2*real_part(cos_integral(d*x
+ a*d/b))*tan(1/2*a*d/b) + 4*a^2*b^2*real_part(cos_integral(-d*x - a*d/b))*tan(1/2*a*d/b) + 2*a^3*b*d - 2*a^2
*b^2*imag_part(cos_integral(d*x + a*d/b)) + 2*a^2*b^2*imag_part(cos_integral(-d*x - a*d/b)) - 4*a^2*b^2*sin_in
tegral((b*d*x + a*d)/b) - 12*a^2*b^2*tan(1/2*d*x) - 12*a^2*b^2*tan(1/2*c))/(b^7*x^2*tan(1/2*d*x)^2*tan(1/2*c)^
2*tan(1/2*a*d/b)^2 + 2*a*b^6*x*tan(1/2*d*x)^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + b^7*x^2*tan(1/2*d*x)^2*tan(1/2*c
)^2 + b^7*x^2*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + b^7*x^2*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + a^2*b^5*tan(1/2*d*x)^2
*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + 2*a*b^6*x*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*a*b^6*x*tan(1/2*d*x)^2*tan(1/2*a*d/
b)^2 + 2*a*b^6*x*tan(1/2*c)^2*tan(1/2*a*d/b)^2 + b^7*x^2*tan(1/2*d*x)^2 + b^7*x^2*tan(1/2*c)^2 + a^2*b^5*tan(1
/2*d*x)^2*tan(1/2*c)^2 + b^7*x^2*tan(1/2*a*d/b)^2 + a^2*b^5*tan(1/2*d*x)^2*tan(1/2*a*d/b)^2 + a^2*b^5*tan(1/2*
c)^2*tan(1/2*a*d/b)^2 + 2*a*b^6*x*tan(1/2*d*x)^2 + 2*a*b^6*x*tan(1/2*c)^2 + 2*a*b^6*x*tan(1/2*a*d/b)^2 + b^7*x
^2 + a^2*b^5*tan(1/2*d*x)^2 + a^2*b^5*tan(1/2*c)^2 + a^2*b^5*tan(1/2*a*d/b)^2 + 2*a*b^6*x + a^2*b^5)

Mupad [F(-1)]

Timed out. \[ \int \frac {x^2 \sin (c+d x)}{(a+b x)^3} \, dx=\int \frac {x^2\,\sin \left (c+d\,x\right )}{{\left (a+b\,x\right )}^3} \,d x \]

[In]

int((x^2*sin(c + d*x))/(a + b*x)^3,x)

[Out]

int((x^2*sin(c + d*x))/(a + b*x)^3, x)